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Abstract. We investigate how the density of baryonic and cold dark matter, the 
density of dark energy and the value of the Hubble parameter at the present time 
influence the propagation of ultrahigh energy protons in the nearby Universe. We take 
into account energy losses in the cosmic microwave radiation, the only one relevant for 
protons above 10 18 eV, and we explore the dependence of Greisen-Zatsepin-Kuz'min 
(GZK) horizon on the cosmology. We investigate several cosmological scenarios, 
from matter dominated to energy dominated ones, and we consider the impact of 
uncertainties in the Hubble parameter in a A— Cold Dark Matter (CDM) Universe, 
estimated from recent observations, on the GZK horizon. The impact of the (unknown) 
extragalactic magnetic field on our study is discussed, as well as possible probes of the 
Hubble parameter attainable by current and future experiments. 
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1. Introduction 

The nucleonic component of extragalactic ultrahigh energy cosmic rays (UHECRs) 
above 100 EeV (1 EeV = 10 18 eV) could be subjected to a strong attenuation because 
of the cosmic microwave background radiation, as it was first noted by Greisen [I] 
and, independently, by Zatsepin and Kuz'min [2]. The main consequence of such a 
predicted Greisen-Zatsepin-Kuzmin (GZK) effect should be to constrain the maximum 
propagation distance of nuclei from extragalactic sources. Moreover, if UHECRs above 
~ 50 EeV are mainly produced by extragalactic sources distributed on cosmological 
distance, then a flux suppression is expected at the highest energies. Recently, the Pierre 
Auger Collaboration [3] and the HiRes Collaboration [4j reported the experimental 
evidence of the suppression of the UHECRs spectrum with a statistical significance 
of about six and five standard deviations, respectively. Although such results do not 
provide a definitive evidence for the existence of the GZK effect, because it can be related 
to a change in the shape of the injection spectrum at the sources, they agree with what 
it is expected from the GZK effect for protons or iron nuclei. More recent observations 
reveal a suppression of the spectrum above 40 EeV with significance greater than 20 
standard deviations [5]. 

Within this study we investigate the impact of cosmology on the propagation of 
protons, and, in particular, on their GZK horizon, i.e. the distance within which 90% 
of observed protons above a certain energy threshold are expected to be produced. In 
fact, protons with energy above 1 EeV lose energy because of photopion production due 
to baryonic resonances, strange particle and multipion production, and pair production 
by interacting with ambient photons of the cosmic microwave background, the infrared 
and optical backgrounds participating marginally [U [2j El [7J [91 [101 [11]. In the case of 
UHE heavier nuclei the existence of additional energy-loss processes in the microwave, 
infrared and optical backgrounds, due to the giant dipole resonance, the quasi-deuteron 
effect and the photofragmentation, drastically reduce the corresponding propagation 
distance, even if by an accident of nature the GZK horizon for both iron nuclei and 
protons is approximately the same [8]. In the following, we will only consider protons. 
In this case, the continuous energy-loss approximation can be safely adopted for the 
study of the GZK horizon [8] (and Refs. therein). 

However, it is worth remarking that many sources of uncertainty in the GZK 
horizon of protons have been already pointed out. First, a genuine limitation to the 
the study of the GZK horizon signatures is provided by the relative uncertainty in the 
energy of UHECR events due to the experimental resolution. A second well-known (but 
removable) limitation is given by the continuous energy loss approximation adopted by 
some authors to simplify calculations. A difference of about 10% in the estimation of the 
GZK horizon of protons emerges if more accurate Monte Carlo simulations are carried 
out instead of a simplified analytic treatment |12] . Within our study, we will show that 
the uncertainty in the values of cosmological parameters, as the Hubble constant at the 
present time, introduces a further independent uncertainty in the horizon. 
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Although several models for production mechanisms of UHECRs are available 
P31 Q3] (and Ref. therein), [151 HSl HZl HH1 HH1 [20l [2U [22] , it is generally accepted that 
the candidate sources are extragalactic and trace the distribution of luminous matter 
on large scales [23J. In particular, it has been shown that correlation with possible 
high redshift sources is unlikely [24J, whereas compact sources are favored [251 26J: 
the recent result reported by the Pierre Auger Collaboration, from observations in 
the southern hemisphere, experimentally supports the latter claim, showing an high 
correlation between the observed data and the distribution of nearby active galactic 
nuclei (AGN) [27J [28j [29]. However, the result has not been confirmed by the HiRes 
Collaboration, from observations in the northern hemisphere [30] . 

Particles produced by a single source S with injection energy E^ are subjected to 
energy losses, and they are detected on the Earth with a degraded energy Ef < Ei. In 
any study devoted to explain the current anisotropy signal or the amount of correlation 
with a particular catalog of sources, it is of fundamental importance to take into account 
a model for UHECRs injection and propagation effects. The expected number of 
particles from a given source at a distance z and with injection spectrum Q(E), can be 
reasonably considered to be proportional to its luminosity £, to the factor z~ 2 and to 
the attenuation factor 



accounting for energy losses. In particular, such a function, when properly normalized, 
estimates the probability that particles, produced at redshift z with initial energy greater 
or equal than Ei, might be detected with energy Ef. The surviving fraction of particles 
with energy above Ef which has been produced in the nearby Universe within a certain 
distance z, Qgzk(z; Ef), depends on ugzk(z; Ef) and will be defined further in the text. 
In this study we consider a pure power-law injection spectrum Q(E) oc E~ s , where s is 
the injection index, as in recent studies [3TJ E2J [331 Ell- 
in this work we investigate the influence of cosmology on the function Qgzk(z] Ef). 
We will consider very different models of the Universe, from flat to curved ones, 
and different values of main cosmological parameters. In particular, we will discuss 
the impact of uncertainty in the Hubble parameter at the present time, according 
to the ACDM model of the Universe and to the experimental constraints obtained 
from recent WMAP observations [35J. However, it is worth remarking that our results 
will be obtained under the assumption that i) the distribution of UHECR sources is 
homogeneous and ii) the luminosity of such sources, emitting a proton-only composition, 
is known and equal for all sources. 

In Sec. [2] we briefly present the main cosmological parameters involved in our 
analysis and the values adopted to investigate some representative models of the 
Universe. In Sec. [3] we discuss in detail the energy-loss processes that UHE protons 
are subjected to and we show how such processes, together with cosmological models, 
are taken into account in the definition of Qgzk(z] Ef). We dedicate part of this section 
to estimate the effect of extragalactic magnetic field on our analysis: we show how the 




(1) 
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impact of its turbulent component on the GZK horizon is small and how it can be 
safely neglected in this study, where only the propagation of protons is considered. In 
the case of heavy nuclei, the extragalactic magnetic field has a dramatic influence on 
their propagation, invalidating their possible use in this work as a cosmological probe. 
Finally, we investigate and discuss the influence of cosmological parameters on the GZK 
horizon, playing a fundamental role in the search of sources of UHECRs. 



2. Cosmological scenarios 



We start by considering the Einstein equation to describe the gravitational field 
(including the term containing the cosmological constant A) in the classical General 
Relativity framework. Under the assumptions of an isotropic and homogeneous 
Universe, we consider the Friedmann- Robert son- Walker (FRW) metric with the 
parameter k accounting for the spatial curvature: k — —1 denotes an open metric, 
u = 0a flat metric and k — 1 a closed metric. Indeed, we consider the Universe as a 
perfect fluid: within such assumptions, the Einstein equation leads to the well known 
Friedmann equations. By introducing the critical density, defined as g c = 3H 2 /8nG, 
where H is the time-dependent Hubble parameter and G is the Newton gravitational 
constant, Friedmann equations can be rewritten as a function of a dimensionless density 
parameter, suitable for the comparison of different cosmological models. 

The density parameter accounts for the matter and the energy in the Universe, 
and can be parametrized as the sum of different contributions. In the standard ACDM 
model, there are some contributions to Q: due to baryonic matter, Q c due to cold dark 
matter, JIa due to dark energy, Q r due to radiation and Q K for the spatial curvature. If 
we define the redshift z by 1 + z = a _1 (i), being a(t) the scale factor, the first Friedmann 
equation can be written in terms of z and of density parameters as 

^fi = a r (i + z y + q m (i + z f + n k (i + zf + q a , (2) 

where Q M = Q b + Q c is the total density of matter and H is the Hubble parameter at 
the present time. By taking into account that the radiation density is important only in 
the early Universe, i.e. at high redshifts, whereas in practice it is negligible in the late 
Universe, the constraint Qm + fi K + = 1 for the density parameters can be obtained 
from very general considerations. Finally, changes in the expansion rate of the Universe 
are described by the deceleration parameter 



from which 



q(z = 0) =q = ]p M - Q A (4) 



at the present time. The parameter q Q and the density parameters described above are 
varied to reproduce very different cosmological models. However, it is worth remarking 
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that we are under the assumptions of an isotropic and homogeneous Universe in the 
approximation of perfect fluid. 

Although the most recent observations based, for instance, on WMAP 
measurements |35j indicate that the curvature of the Universe is very close to flat, 
within the present study we are also interested in investigating the impact of non-flat 
cosmology on the GZK horizon of UHE protons. 

In a flat Universe, we consider the deceleration parameter q constrained by 
VLm + = 1, whereas in a curved Universe, we consider the curvature parameter 
Q K constrained by Q K = 1 — Qm — ^a- By varying two among the three parameters, 
we investigate different models of the Universe. In Tab. CD and Tab. [2] the values of the 
parameters are summarized for flat and closed models, respectively, that will be adopted 
in the successive analysis. 



Name £Im Qo Description 



ACDM 0.272 0.728 -0.592 

EdS 1 0.5 

AFU1 1 -1 

AFU2 I | -0.5 



Standard cosmological model 
Einstein-de Sitter model 
(matter-dominated) 
Decelerating flat Universe 
(vacuum energy-dominated) 
Decelerating flat Universe 



Table 1. Models for flat Universes considered in this study and corresponding to 
different values of the density of matter, the density of dark energy and the deceleration 
parameter. 



Name 


Qm 


^A 


f2 K 


Description 


FLOl 


0.27 


0.93 


-0.2 


Friedmann-Lamaitre open 


FL02 




1.23 


-0.5 




FL03 




1.43 


-0.7 




FLC1 


0.27 


0.53 


0.2 


Friedmann-Lamaitre closed 


FLC2 




0.23 


0.5 




FLC3 




0.03 


0.7 





Table 2. Models for curved Universes. 



3. Propagation of UHE protons 

The propagation of cosmic rays is generally treated as a diffusive process described by 
the Ginzburg-Syrovatskii transport equation [36J, and later reviews. Such equation is 
rather complicated, allowing an analytical solution in few cases under very restrictive 
assumptions. In the case of UHECRs propagating in the extragalactic space, the 
transport equation can be simplified by assuming time-independent diffusion coefficient 
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and energy losses, defining a static Universe [37J ESI ESI ED] • A more general approach, 
including such time dependence for ultrarelativistic particles diffusing in an expanding 
Universe from a single source, has been recently proposed [41] . Other important 
ingredients required to correctly describe the propagation of UHECRs are cosmic 
background radiations and magnetic fields. As we will see in the following, background 
radiations are responsible for the energy loss of UHECRs, whereas magnetic fields 
influence their trajectory and may have a non-negligible impact on the propagation 
time and on the probability to reach the Earth. For the purpose of our study, we focus 
on the term 

^ = -[H(t)E + b int (E,t)}, (5) 

describing energy-loss due to the expansion of the Universe and to proton interactions 
as a function of time [JT]. In the following, we will consider the energy- loss rate as a 
function of the redshift z, defined by 

_ x 1 dE ^, dt 
E dz ' dz' 



in a Friedmann Universe, where 
dt 1 



[n M (i + z ) 3 + n A + n K (i + z) 2 \ 



dz H (l + z) 

is the factor accounting for the cosmological expansion [121 H21 E]. The function 
(3(z, E) is related to the cooling rate of protons and it depends on the particular 
energy- loss process considered. In the energy interval of our interest, three main energy- 
loss processes can be considered: i) adiabatic, due to cosmological expansion, ii) pair 
production and iii) photoproduction caused by the interaction with cosmic microwave 
background (CMB) radiation, also known as the Greisen-Zatsepin-Kuzmin (GZK) effect 
[H [2] . The energy loss due to inverse Compton effect is negligible, for energies greater 
than 10 17 eV, and it will not be considered in the following. 

We adopt the standard black body model with temperature T ~ 2.725 K for the 
CMB. If e denotes the photon energy in the observer's rest frame, the photon energy 
density at the present time is defined by 

x- 1 



n(e) e 2 

' exp 



e 2 7T 2 (hc) 



k B T 



The adiabatic term, accounting for energy loss rate due to the expansion of the Universe, 
is given by 

(3 rsh {z) = H p M (l + zf + Ov + tt K (l + zf] 1 , (6) 

obtained from Eq. §2§. Parameterizations for pair and photomeson production are 
discussed in the following, whereas the impact of magnetic fields on our study will 
be discussed at the end of this section. 
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3.1. Pair production 

In the rest frame of the proton, pair production process occurs at the threshold energy 
2m e c 2 ~ 1 MeV and it plays an important role only when CMB is considered, the CIB 
participating marginally [6]. We treat the process as a continuous energy loss, because 
the loss per interaction is very small, about 10~ 3 x E. We consider the energy loss rate 
j3 e ± (E) accounting for the pair production, because of the Bethe- Heftier interaction with 
ambient photons with density n(e), defined jl5] by 

'm e c 2 \ ip(£) 



(3 e± (E) oc ar 2 e Z 2 (m e c 2 ) 2 jT \n 



2 



where 7 ~ E/m p c 2 is the Lorentz factor of the proton, m e is the electron mass, a = e 2 / he 
is the fine-structure constant and r e = e 2 /m e c 2 is the classical electron radius. In the 
case of a black body radiation with temperature To, as in the case of CMB, the energy 
loss obtained in the Born approximation becomes 



ar 2 e (m e c 2 k B T ) 2 c 
M*0 = -^3-^ /(„), 



with 



roc 

f\v) = v 2 J dt<p(£) (e«« - 1) 



-1 m e c 2 



2jk B T 

where k B is the Boltzmann constant. Higher order terms of the Born approximation, 
proportional to (Zav±/c) m , where m is the number of interactions with the Coulomb 
field, should be taken into account in the case of nuclei heavier than protons. In fact, for 
Z > 1 the symmetry between produced electron and positron breaks down. Blumenthal 
suggested to correct the rate through the Sommerfeld factor [15], although it is only 
valid in the non-relavistic limit. For nuclei with Z > 1, this correction does not agree 
with experimental data and a better correction is required |46j . By taking into account 
the evolution factor, after some algebra, we obtain 

/M^#rv r r® , (7) 



E 3 



exp 



1 



where the auxiliary function </?(£) is parametrized in Ref. [17], masses are in units of 
eV/c 2 and 



■h 2 2 

A e ± = - I J w 3.44 x 10- 18 EeV 3 s-\ 
= mjn^ ^ E y 
2k B T 

3.2. Photonuclear interactions 

The probability of UHE protons to interact with CMB photons rapidly increases with 
proton energy. In fact, in the proton's rest frame the interaction is equivalent to a 
collision with a high energy photon with energy e'. When the energy e' equals at least 
the pion mass m^c 2 m 140 MeV, the proton undergoes photomeson production and loses 
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energy. Such a process is known as Greisen-Zatsepin-Kuzmin effect and dominates above 
50 — 60 EeV [H The two main channels for the interaction, close to the threshold 
energy, are 

p + 7 — ► A(1232MeV) — > n + n + , 

n — > p + e" + v e . 

involving the resonance A(1232MeK). At highest energies, heavier resonances and 
multipion production channels are likely. The energy-loss rate due to the interaction on 
CMB, the dominant one above 50 EeV, is given by 

,2 poo „/ \ r 2e- E 



k B T m 2 p 
2vr 2 n £ 2 



deK.(e)a(e)e x In 



1 — exp 



2Ek B T 



where m p is the proton mass in units of eV/c 2 , a(e) is the cross-section for pion 
production in terms of the photon energy e and /C(e) is the inelasticity factor [HJ 119] . 

Just above the threshold, baryonic resonances dominate and protons are subjected 
to photomeson production, mainly through the A(1232)-baryon resonance, whereas 
heavier resonances (up to A(1950)-baryon) play a more marginal role. The cross-section 
for baryonic resonances is parametrized by 

4 

oWe) = ^ViCr L (e;ei,rj) 

i=l 

where o~l is the Lorentzian function, (e, (GeV), Tj (GeV), (fib)) = (0.34,0.17,351), 
(0.75,0.50,159), (1.00,0.60,21) and (1.50,0.80,26) for i = 1,2,3 and 4, respectively. 
For all other processes participating in photomeson production, including multipions 
(MP) or direct particle production involving ir, rj, A, p, oj and strange-particle channels 
(RP), we use Rachen's parameterizations [46J. In Fig. [1] are shown the cross-sections for 
the discussed pj interactions, separately, and the total cross-section. 

By taking into account the evolution factor, the energy-loss rate for photomeson 
production is given by ^{z,E) = (1 + zf^(z = 0,E(1 + z)) for CMB. Finally, the 
total energy-loss rate is defined as 

E~^ = ~Tz X ^ rsh ^ + + ^^ E ^ ' ( 8 ) 

taking into account the contributions of the processes discussed in this section. Efficient 
analytical parameterizations of the cooling rates can be considered for the numerical 
estimation of the energy loss. First, we have verified that the recent parameterization 



A w (l + z) s exp 



E < E(z) 



(3 n (z,E) ~ <^ ' "> Y^ E \ ~ y> (9) 

\a(l + z) 3 E>E(z) 

provides an excellent approximation for the energy-loss rate due to photomeson 
production. Here, the function E(z) = 6.86e _a8072 x 10 20 eV ensures the continuity 
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Photon energy, e' (GeV) 



Figure 1. Cross-section for the photomeson production of a proton in the CMB 
radiation. Contributions from baryonic resonances (BR), Rachen's parameterizations 
for direct particle production (RP) and multipions (MP) are shown. Energy has to be 
considered in the rest frame of the proton. 




Figure 2. Left panel: Energy loss length in the CMB, for each process separately and 
for all processes together, in the case of a proton; triangles indicate the values recently 
reported by Stanev |51j . Right panel: GZK horizon as a function of the propagated 
proton energy (triangles) and values obtained by Harari et al [8] (dashed line). 



of /3 n (z,E) and {A n ,B^,C^} = {3.66 x 10-V" 1 , 2.87 x 10 20 eU,2.42 x lCrV -1 } are 
taken from Ref. [50] . 

In Fig. [2] (left panel) is shown the energy-loss length in the CMB, for each process 
separately and for all processes together, in the case of a proton propagating in a ACDM 
Universe. In the energy interval between 1 EeV and 50 — 60 EeV the main energy-loss 
process is the pair production, whereas photomeson production dominates up to the 
highest energies. We show, for comparison, the excellent agreement with the total loss 
length recently reported by Stanev, for a proton propagating in the CMB [51] . 

The solution of Eq. (jSJ) is the degraded energy Ef of the proton with initial energy 
Ei, after the propagation from the source at redshift z to the Earth (z — 0). Under the 
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assumption of a power-law injection spectrum, we have the surviving probability 

s 
E 



s — i r 

uj G zk(z; Ef) = / E~ s dE, (10) 

/ JEiUEf) 



while assuming equal intrinsic luminosity and homogenous distribution of sources (i.e. 
their number at redshift z is proportional to z 2 ), we obtain the surviving flux defined 
by 

, N rdz'f£ (z , Ef) E-*dE 
Q G zk{z; E f ) - -00-77700 j?- s jj?> C 11 ) 

Jo dz J Ei (z';E f ) E dE 

for the probability of detecting at z = a proton emitted with E > Ef farther than 
a distance z. The energy Ei (z; Ef) of the injected proton is estimated by evolving 
Eq. OH]) backward in time [HI [33] . The GZK horizon R is the distance such that 
1 — &>gzk{R, Ef) = 0.9. In Fig. [5] (right panel) is shown the GZK horizon as a function 
of the propagated proton energy and we show, for comparison, the good agreement with 
results obtained by Harari et al [8]. 



3. 3. Extragalactic magnetic field 

It is worthwhile discussing in this section the impact of the extragalactic magnetic field 
(EMF) on the probability Qgzk(z] Ef). Unfortunately, a direct or indirect measurement 
of the EMF is still missing, and only bounds to its r.m.s. strength B rms have been 
estimated. However, bounds on the EMF strength also depend on the field correlation 
length £, which is also unknown. 

Upper limits of the order of 10 -9 G (1 nG) on the intensity of the EMF have been 
measured through i) the Faraday rotation in the polarized radio emission from distant 
quasars [52] and ii) the characteristic distortions that it induced on the spectrum and 
the polarization of CMB radiation [53]. More recent analyses, based on the recent 
WMAP measurements [35], provide the more stringent upper bound of ~2 nG on the 
present value of the cosmic magnetic field of primordial origin, more than one order of 
magnitude smaller than previous estimates 

Intergalactic magnetic field might also be structured inside and around clusters or 
groups of galaxies, with filaments extending over few Mpc, as shown, for instance, in 
recent detailed simulations [55]. The topology of such a structured magnetic field would 
have a non-negligible impact on the trajectories of UHE protons and, of course, on the 
average deflections they experience in the case of EMF with regular structures above 
200 kpc [56]. Additionally, longitude-averaged X-ray emission observed with ROSAT 
near 0.65 keV and 0.85 keV towards the center of the Galaxy, are in agreement with 
a Galactic wind thermally-driven by cosmic rays and hot gas [57, 58J. In our study, 
we can neglect the effect of such a magnetic wind, because it is expected to have a 
non-negligible impact for protons below 60 EeV or heavier nuclei. Moreover, we neglect 
the case of a structured EMF because of the lack of direct or indirect measurements 
about its structure. 
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Figure 3. Relative difference between GZK horizons estimated with and without 
EMF. Two values of the r.m.s. strength of the EMF are considered, namely 1 
and 10 nG, as well as two different injection index, namely 2.4 (left panel) and 2.7 
(right panel). Error bars indicates only the statistical uncertainty on SR/R. Lines 
corresponding to ±5% bands and SR/R = arc shown for reference. 



Here, we consider only the case of a turbulent EMF characterized by B rms and I. 
In this case, protons propagate randomly following a brownian trajectory. The average 
deflection they are subjected to can be parameterized by 

e (E , D) * o.8» (^) " yy ' (^y 4 (^) , w 

in absence of energy losses [14J, with D zc/Hq for distances considered in the present 
study. Such a deflection implies an average time delay 

10 20 ev) (lO Mpc) (lO- 9 g) kyr ^ 

relative to the rectilinear propagation with the speed of light. Let us consider a 
propagation distance equal to the GZK radius, i.e. D = R: the difference between 
the trajectory corresponding to the propagation in absence of magnetic field and the 
brownian trajectory of the proton propagating in the EMF is 5R cr, leading to 



<W> _ E \~ 2 ( R W £ \(B T x 2 



i?-°- 5xl0_ [l^Vj { 10Mp~c ) \TMp~c ) [lO^G ) ' (M) 
where we have omitted to explicitly report the dependence of R and 6R on the energy 
E, for simplicity. We estimate the bound to B rms \/i in order to have a negligible impact 
on the GZK horizon, i.e. 5R/R < 5%, by 



I B rms 3 / E \( R \- 

1 Mpc 10- 9 G < 2 VlO 20 eVJ \10 Mpc) ' [ } 

For a UHE proton with E = 60 EeV the GZK horizon is R w 180 Mpc (for 
Ho = 70.4 km/s _1 /Mpc): if we consider the propagation in an EMF characterized by 
i = 1 Mpc, any r.m.s. strength B rms < 4.5 nG will not significantly affect the horizon. 
For a proton with higher energy as 95 EeV, such a bound reads B rms < 11.5 nG, 
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i.e. even more intense EMFs have a negligible impact on the result. For instance, the 
estimated upper bounds are still conservative even if larger correlation length as 16 Mpc 
are considered, reducing to ~1 nG and ~3 nG, respectively. 

We have investigated the validity of such an argument by performing several Monte 
Carlo realizations of protons propagating in the Universe, with and without EMF. For 
such a purpose, the well known propagation software CRPropa [59J has been adopted: 
two different values of the r.m.s. strength of the EMF, namely 1 and 10 nG, as well as 
two different values of the injection index, namely 2.4 and 2.7, has been considered for 
this study. In any case, the correlation length has been fixed to I = 1 Mpc. The main 
advantage in using CRPropa is that relevant energy-loss processes are considered during 
the propagation in a turbulent magnetic field: in this case we can quantify the average 
relative deviation 5R/R in a more realistic scenario with respect to Eq. (I14p . where 
energy loss has been switched off to estimate 9(E,D) and t(E,D). In Fig. [3] is shown 
the result of such a study. It is evident that the the injection index has a negligible 
impact on the relative deviation 5R/R. For both B rms ~ 1 nG and B rms ~ 10 nG, the 
resulting relative deviation is within 5% at any energy, in agreement with expectation. It 
is worthwhile noticing that in any case, the value of 5R/R is non-positive, i.e. the GZK 
horizon for the propagation in an EMF is smaller than the horizon for the propagation 
without magnetic field, as axpected. Moreover, it is relevant to remark that for the 
weak EMF (1 nG) 5R/R « 0, whereas for the strong EMF (10 nG) SR/R approaches 
zero above 90 EeV: the EMF has no impact on the GZK horizon of UHE protons. 

From such a result, we can also deduce that above 60 EeV the effect of energy 
loss on the propagation in an EMF is negligible: hence, the parameterizations given by 
Eq. ( 112"]) and (TIB"]), although they have been obtained by considering no energy losses, 
can be safely used in the scenarios previously described. 

Such a result is rather conservative: in fact, the most up-to-date bounds on the 
correlation length and the r.m.s. strength of the EMF are much lower of the values 
considered in this study [60] . 

4. Analysis and discussion 

In Fig. H] is shown the function Qgzk{z; Ef) for some values of the energy threshold 
Ef = Ethri as a function of the distance. By assuming the most recent values of ACDM 
model parameters {Q b = 0.0456, Q c = 0.227, Q A = 0.728 and H = 70.4 km/s^Mpc) 
[35], the value of the injection index is fixed (s = 2.7) while the energy threshold E t hr 
is varied (left panel): as expected, the GZK horizon decreases by increasing the energy 
threshold. In the right panel of Fig.H] is shown the same function for two fixed values 
of the energy threshold, E t h r = 50 EeV and E t h T — 100 EeV, while the injection index 
is varied: although differences among curves corresponding to different values of the 
injection index are not exaggerated as in the previous case, we find that they increase 
by decreasing the energy threshold E thr . Although observations suggest an injection 
index between 2.2 and 2.6, depending on the underlying assumptions, we have extended 
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our study to a broader range, namely from 2.0 to 2.7, because of interest for recent 
studies [MUSHES]. 

Cosmological parameters as the density of matter, the density of dark energy and 
the Hubble parameter at the present time, should influence the function Ei(z;Ef), 
discussed at the end of the previous section, and by consequence the probability 
Qgzk{z; Ef). From Eq. (jSJ it is clear that the factor involving such parameters is only 
—dt/dz. 

For small values of z, corresponding to the nearby Universe, it is not expected a 
significant difference among cosmological models. At the highest energies, above 50 
EeV, we have found (see Fig. [2j left panel) that photomeson production dominates 
energy losses, thus only the term —f3 w (z,E) x dt/dz is important. By considering 
two extreme scenarios, as matter-dominated (Qm — 1; = 0) or energy-dominated 
(Qm = 0; = 1) Universes, the energy-loss term is modified at most by 10%. 

In Fig.|5]is shown the function Qgzk(z; Ef) as a function of the redshift, for both 
flat and curved models listed in Tab. [Hand [2], respectively. With no regards for the value 
of the injection index s, functions corresponding to very different cosmological models 
do not differ significantly below the GZK horizon, which is not influenced by the choice 
of a particular cosmology. The difference becomes significant toward higher redshifts 
where, however, the surviving probability for a proton is very small. 

Indeed, in Fig. [6] is shown the same function by assuming ACDM model values for 
density parameters [35], protons with energy threshold E t hr = 60 EeV and different 
values of the Hubble constant H at the present time. For simplicity, in the following 
we will not explicitly specify the unit of Ho, that should be considered km/s _1 /Mpc. 
Left panel shows the case for s = 2.0, whereas right panel shows the case for s = 2.7. 
It is evident that the value of the Hubble parameter significantly affects the function 




Figure 4. Surviving flux, denned by Eq. ([TTj) . due to cosmological effects described in 
the text, in the case of a proton and assuming ACDM model parameters. Left panel: 
the value of the injection index is fixed (s = 2.7) while the energy threshold Ef is 
varied; Right panel: the value of the energy threshold Ef is fixed (two cases, E = 50 
EeV and E — 100 EeV) while the injection index is varied. 
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Redshift, z Redshift, z 

Figure 5. Surviving flux, denned by Eq. (fTTj) . due to cosmological effects described in 
the text, in the case of a proton with E t hr = 60 EeV and assuming different models of 
the Universe, for two values of the injection index (s = 2.0 and s = 2.7). Left panel: 
flat models described in Tab.Q] are considered; Right panel: curved models described 
in Tab. [2] are considered. 



^gzk(z; Ef), as already suggested by previous studies [49j [31]. If the 3a uncertainty 
in the value (Hq = 70.4) obtained from the ACDM model is taken into account [35], as 
shown in both panels of Fig-EJ we find that differences between curves corresponding to 
70.4 + 3(7 and 70.4 — 3a are negligible (« 5%) only for z < 0.018, as shown in the inset 
of the right panel, for s = 2.7. From the lower bound curve (Hq = 70.4 — 3er ~ 66) 
the estimated GZK horizon is z ~ 0.039 (~ 156 Mpc), whereas for the upper bound 
curve (H£ = 70.4 + 3a « 75) it is z « 0.044 (« 198 Mpc), with a relative difference 
of about 11% in redshift and about 21% in distance. A relative difference in redshift of 
about 100% between the corresponding surviving functions is reached around z = 0.037, 
tending to increase with the redshift. 

In Fig. [7] are shown the relative differences between curves corresponding to 
surviving functions for Hq and Hq , by varying the energy threshold and the injection 
index, s = 2.0 (left panel) and s = 2.7 (right panel). By increasing the energy threshold 
the relative difference increases, although it keeps smaller than 15% below the GZK 
horizon. All significant differences are found to be above the GZK horizon, beyond 
which only 10% of protons output by their sources are able to reach the Earth. In the 
following we will discuss two simple applications of the arguments discussed so far, in 
order to show how current and future experiments could probe the value of the Hubble 
parameter. 

4-1. Application #1: energy spectrum 

As a first practical application, we consider the propagation of 10 s UHE protons 
produced with CRPropa [59] by sources homogeneously distributed within 2 Gpc. It 
is worth remarking that CRPropa performs propagation of particles with respect to 
distance, for a specified astrophysical and cosmological scenario. Hence, in this section, 
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Figure 6. Surviving flux, as in Fig. [31 assuming ACDM model values for density 
parameters and protons with energy threshold E t hr = 60 EeV, for two values of the 
injection index (s = 2.0 and s = 2.7). Different values of the Hubble parameter Ho (in 
km/s _1 /Mpc units) at the present time are considered. The inset in the right panel 
shows the relative difference between curves corresponding to Hq = 70.4 + 3a s» 75 
and Hq = 70.4 — 3a ps 66, as a function of the distance. 




Figure 7. Relative difference between curves corresponding to probability functions 
for Hq and Hq in Fig. [5J for different values of energy threshold in the case of s = 2.0 
(left panel) and s — 2.7 (right panel). 



we will report the results with respect to distance. Moreover, we consider neutrinos 
produced by UHE protons and propagated by CRPropa. The injection spectrum at 
source is considered to be oc E~ 2 ' 7 , with a maximum energy of E max = 1000 EeV. The 
propagation of protons that have reached an energy smaller than 10 EeV, because of 
the interactions with CMB photons, is not followed. In particular, we consider three 
astrophysical scenarios, with different values of the Hubble parameter at the present 
time, namely H = 70.4, Hq and Hq . In Fig. [8] (left panel) is shown the number of 
entries versus the distance at which GZK neutrinos are generated, together with the 
flux of GZK neutrinos at Earth, for each astrophysical scenario (right panel). While for 
distances below ~ 1 Gpc, differences in the three scenarios are small, at large distances, 
comparable to those of the most far sources, a larger number of neutrinos is created for 
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increasing values of the Hubble parameter. 

The interaction rate is expected to decrease for increasing values of Hq, because 
of its dependence on the cosmological factor dt/dz, which in turns depends on the 
inverse of Hq. Although differences are small, such an expectation is confirmed at the 
lower distance (see left panel of Fig.[8]). On the other hand, the number density of 
background photons increases with distance, as well as their energy: such a behavior 
tends to dominate the previous one for distance above ~ 1 Gpc. The differences in 
the distance at which GZK neutrinos are produced, have a direct impact on the energy 
spectrum of protons (<& P (E)) and neutrinos (<& U (E)) observed at Earth. In the case of 
neutrinos, Q V (E) is shown in the right panel of Fig.|HJ 

The relative differences between the flux of protons and neutrinos corresponding to 
Hq and Hq, are shown in Fig. [91 The effect of varying the Hubble parameter is a simple 
diagonal shift in the flux. Let 5<&(E)/<&(E) indicate the relative difference of the energy 
spectrum corresponding to scenarios with H with respect to scenario with Hq. In the 
case of protons, where a simple power-law flux is expected up to 40 — 50 EeV, such a shift 
induced by Hq should produce a rather constant 8<& P (E) /<& P (E) versus energy. Such an 
expectation is confirmed by the result shown in the left panel of Fig.[9j where constant 
differences of the order of 4-5% are found at any energy below ~60 EeV. We argue 
that such differences can be quantitatively explained by the different energy-loss rates 
corresponding to the values of Hq and Hq. In fact, the energy- loss rate is proportional 
to the cosmological factor —dt/dz, and, by consequence, to the inverse of the Hubble 
parameter at the present time. It can be shown that the diffuse spectrum reflects such a 
proportionality to the inverse of H because of its direct dependence on the cosmological 
factor (see, for instance, Ref. |63j). Hence, the expected difference in the flux of UHE 
protons approximately reduces to (H — Hq)/ Hq, providing a ~ 5% alteration of the 
flux in the case of Hq and a ~ —5% alteration of the flux in the case of Hq, in good 
agreement with our finding. 

In the case of neutrinos, the more complicated dependence on energy of Q v is 
expected to produce non-constant relative differences 6<& U (E) / <& U (E). The results shown 
in the right panel of Fig. [9] confirms such an expectation: relative differences range from 
0% to pa 10%. Such a result suggests that experimental evidences of this deviation should 
be explored below 10 16 eV, around 10 17 eV or above 10 19 eV. However, it is worth 
remarking that a method to exploit such a deviation to probe the Hubble parameter is 
still under investigation and it will be the subject of successive studies. 

4-2. Application #2: clustering 

Motivated by the recent correlation between the arrival directions of UHECRs detected 
with the Pierre Auger Observatory and AGN [29J , as a second practical application, we 
consider the distribution of sources corresponding to the position of AGN in the nearby 
Universe (up to z = 0.047), reported in the SWIFT-BAT 58-months catalog |64j. In 
particular, we consider two scenarios: i) we assume equal intrinsic luminosity for all 
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Figure 8. UHE protons produced by sources homogeneously distributed within 2 Gpc 
(see the text for further detail); injection spectrum follows a power-law with spectral 
index s = 2.7 and E max — 1000 EeV. Left panel: number of entries versus the distance 
at which GZK neutrinos are generated, by varying the Hubble parameter at present 
time; Hq = 70.4, Hq and Hq are considered. Right panel: corresponding fluxes of 
GZK neutrinos at Earth. 




Figure 9. Same simulation setup as in Fig. [8] relative difference between fluxes at 
Earth corresponding to Hq and Hq , with respect to Hq, in the case of protons (left 
panel) and neutrinos (right panel). 



AGN and ii) the intrinsic luminosity of each AGN is taken into account. If C indicates 
luminosity of an AGN, and by assuming no source evolution, the probability to get an 
event from such a source is proportional to C z~ 2 ugzk( z , E thr ). Thus, protons are then 
propagated in a ACDM Universe until they reach the Earth. We consider only UHECRs 
with energy above 100 EeV and with arrival direction lying in the field of view of the 
Pierre Auger Observatory, whose non- uniform exposure is taken into account, as well as 
its angular uncertainty of 0.8°. The effect of EMF is also taken into account, smearing 
the direction around the source by sampling a Fisher-von Mises distribution, i.e. the 
Gaussian counterpart on the sphere. The spreading angle is given by Eq. ( 1X21) in the 
case of r.m.s. strength B rms = 2 nG and correlation length £ = 1 Mpc, according to 
the most recent upper bounds [M]- Additionally, according to the result reported by 
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Figure 10. Expected clustering signal, as a function of the angular scale (in deg units), 
from a sky of N = 200 protons with E > 100 EeV (in the field of view of Pierre Auger 
Observatory) and for values of the Hubble parameter considered in Fig. [6] Sources 
of 44% of events are AGN within z = 0.047 in the SWIFT-BAT 58-months catalog, 
whereas the remaining 56% of events are isotropically distributed. Scenarios with 
intrinsic luminosity taken into account (left panel) and not taken into account (right 
panel) are considered. The signal at each angular scale is obtained by averaging over 
10 4 Monte Carlo realizations. 
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Figure 11. Same as in Fig.[T0l for scenarios where intrinsic luminosity of AGN is 
taken into account. The clustering for N = 100 (left panel) and TV = 50 (right panel) 
is considered. 



Pierre Auger Collaboration in the case of the SWIFT-BAT 58-months catalog, the 56% 
of events in the simulated sky are isotropically distributed [29] . 

We investigate the clustering signal averaged over several Monte Carlo realizations 
(10 4 for each astrophysical scenario), by mean of the novel multiscale autocorrelation 
function (MAF) [65], versus the angular scale. Such a function involves the equal-area 
binning of a spherical region of the sky: the number iV of bins defines the angular scale 
of the analysis. Let ipk{®) be the fraction of points from the data set and ipki®) be 
the expected isotropic fraction in the bin B^. The Kullback-Leibler divergence [66J 

a(q) = v KL ty(e)p(e)) = E^(©)iog|^|| (16) 
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quantifies the departure from an isotropic distribution at the scale B. If Ad ata {Q) and 
A iso (Q) refer, respectively, to the data and to an isotropic realization with the same 
number of events, the MAF estimator is defined as the standardized deviation from 
isotropy: 



where (A iso (Q)) and 0^.^(0) are the sample mean and the sample standard deviation, 
respectively, estimated from several isotropic realizations of the data. The MAF 
estimator is not biased against the null hypothesis Ho of an underlying isotropic 
distribution for the data and the estimated chance probability depends on 0. Indeed, 
s(0) follows a half-Gaussian distribution with zero mean and unitary variance, and the 
probability to obtain a maximum value of s(0), at any angular scale G, greater or equal 
than a given value s* is 



providing an analytical expression for the properly penalized chance probability, 
independently on the value of the angular scale G and on the data size. It is worth 
remarking that the angular scale where the chance probability is minimum (i.e. where 
s(G) is maximum) turns to be the most relevant clustering scale and that, for each Q, 
the value of s(G) is an estimation of the amount of clustering at that scale [65J. 

In this study, different values of the parameter Hq are considered, as well as an 
increasing number of events in the sky. The results are shown in Fig. [10] and [TTJ for 
different astrophysical scenarios and angular scales (in deg units). In Fig. [10] we consider 
the cases where intrinsic luminosity of AGN is taken into account (left panel) and not 
taken into account (right panel), for skies of N = 200 protons. In Fig. [11] we focus on 
scenarios where intrinsic luminosity is accounted for, and vary the number of protons, 
to put in evidence the impact of the statistics on the clustering signal. It is evident that, 
for a fixed number of events, the clustering signal increases for increasing values of Hq, 
whereas it decreases for decreasing number of events, as expectecfj]. The dependence 
of the clustering signal on the value of the Hubble parameter at the present time can 
be understood in terms of the surviving probability defined by Eq. ( flOl) and ( [IT]) . In 
fact, we have previously shown that changes in the value of H have a non-negligible 
impact on the surviving flux Qgzk{z)- Such an impact is reflected in the weight function 
^gzk(z) adopted in our simulations for the probability to get an event from a source. 
For a fixed distance, the probability to reach the Earth for UHE protons propagating in 
a ACDM Universe increases for increasing values of Hq, as clearly deducible from Fig. [6] 
Such a behavior favors the clustering around nearby sources, increasing the signal at 
any angular scale. 

A direct comparison, as a function of the angular scale, between the clustering signal 
s(G) obtained from the data and that one obtained from simulations, for different values 



s(9) 



A data (G) - (A lso (®))\ 
^ so (G) 



(17) 





| The statistical power of the method increases with the number of events. 
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Figure 12. Same scenario as in Fig. [IT] with N = 100 events above E t hr = 100 EeV. 
It is shown the cumulative distribution of the clustering function s(G) in the case with 
= 4° (left panel) and = 10° (right panel), for different values of the Hubble 
parameter. For each scenario, 10 4 Monte Carlo realizations have been considered. As 
an illustrative example, we assume that the clustering coefficient obtained from the 
data is 1.4 (solid vertical arrow): in both panels, we show a test with 90% confidence 
level (dotted line) against such a value. See the text for further details. 



of H 0} represents the principal tool for probing the Hubble parameter with our clustering 
approach. The procedure is outlined in the following. 

Let us assume that observations provide a number N of events above a certain 
energy threshold E t h r - We consider an astrophysical scenario with a certain value of the 
Hubble parameter H, by using a catalog of real candidate sources producing a fraction f p 
of UHECRs. The remaining fraction 1 — f p of particles is distributed isotropically in the 
sky. Hence, we simulate a large number of skies with N events in the field of view of the 
observatory (taken into account the non- uniform exposure, if required), and we estimate 
the multiscale clustering s(O) from both the data and the simulations. For a fixed 
angular scale 6o, the statistical test can be performed by comparing the distribution 
of s S i m (Qo] H), obtained from the model, against the value s fc s (6o), obtained from the 
data. For instance, either one- or two-tailed test can be adopted to accept or reject null 
hypothesis that the observation is compatible with the expectation. In order to gain full 
advantage from information provided by the multiscale approach, such a procedure can 
be repeated for any value of the angular scale, and the results can be combined with 
the Fisher's method, its extensions, or any other statistical method. 

As an illustrative example of our procedure, in Fig.[12]we show the case of a one- 
tailed statistical test with 90% CL (indicated by the dotted horizontal line). In absence 
of values of s b s obtained from real data, we assume, for instance, that s Q b s = 1.4 
(indicated by the solid vertical arrow in the figure). Moreover, we consider the same 
astrophysical scenario adopted to obtain the results shown in Fig.dH with N = 100 
events above Ethr — 100 EeV. In particular, we show the cumulative distributions of 
s S j m (O ; H) for three different values of the Hubble parameter, namely 50, 70.4 and 
100 km/s -1 /Mpc, and for two different values of the angular scale, namely 4° and 10°. 
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For such a test, the null hypothesis can not be rejected if the intersection between 
the line corresponding to s b s and the curve corresponding to a scenario with a certain 
Hubble parameter, lies below the CL line. 

Hence, in both examples, it is evident that the null hypothesis should be rejected 
in the cases corresponding to H = 70.4 and H = 100, whereas it can not be rejected for 
H = 50 at 90% CL. 

Hence, such a study suggests that the observation of 100 protons with energy 
above 100 EeV should be sufficient to probe the Hubble parameter. The large statistics 
required for this investigation can be attained in few years of activity by future arrays 
much larger than the Pierre Auger Observatory and with larger exposure, as, for 
instance, JEM-EUSO. However, it is worth remarking that we have considered only 
the most conservative case, with events above 100 EeV. In Fig.[3] is shown that an 
energy threshold of 60 EeV is sufficient enough to observe significant differences in the 
flux, and, indirectly, on the clustering of protons. At such energy threshold, we already 
expect that current experiments as the Pierre Auger Observatory will collect more than 
200 UHECRs within a few years, providing the statistics required for this probe. 

However, it is worth remarking that such a result depends on the assumptions about 
the sources of UHECRs and, by consequence, about their distribution in the nearby 
Universe. In fact, the clustering signal measured in the arrival direction distribution 
of UHECRs is sensitive to the intrinsic clustering of sources and their number density 
|34j . Moreover, the (still unknown) composition of UHECRs plays a significant role in 
the formation of clusters of particles. In fact, both extragalactic and galactic magnetic 
fields have a negligible impact on the deflection of UHE protons. Conversely, magnetic 
fields are expected to bend significantly the trajectories of UHE heavier nuclei, diluting 
the clustering signal at the smallest angular scales and altering the signal at the largest 
ones. In any case, our procedure is robust against the ignorance about the intervening 
magnetic fields and the composition of UHECRs, if the fraction of protons (adopted as 
input to the simulations) is estimated from the observations, as for instance in Ref . [29J , 
and if the remaining fraction of UHECRs is prevalently composed by heavy nuclei. 
Hence, it is clear that without a definitive knowledge of the UHECR source population, 
the inferred value of the Hubble parameter at the present time would depend on the 
underlying assumptions about the candidate sources adopted for the study. 

5. Conclusion 

The GZK effect plays a fundamental role in the search of sources of UHECRs. Within 
the present work we have investigated the influence of cosmology on the GZK horizon 
of extragalactic UHE protons, with energy ranging from 50 to 100 EeV. By considering 
very different models of the Universe, from flat to curved ones, we have shown that 
significant differences among cosmological models appear to be important above the 
GZK horizon, where the surviving probability for the protons is very small. Moreover, 
we have investigated the impact of uncertainty in the Hubble parameter at the present 
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time, according to the ACDM model of the Universe and to the experimental constraints 
obtained from recent WMAP observations. Our results suggest the existence of non- 
negligible differences between the estimated values of the GZK horizon in Universes 
with Hubble parameter H = 70.4 + 3a and H = 70.4 — 3cr, respectively. However, our 
numerical results show that such differences should have a small impact on studies 
involving distances below 250 Mpc, as for instance the recent correlation analyses 
between observed data and the distribution of nearby active galactic nuclei (AGN) 
reported by the Pierre Auger and the HiRes collaborations [271 ES E9[ [30]. Finally, we 
have shown that current and future experiments could probe the value of the Hubble 
parameter by measuring the flux of GZK neutrinos at Earth or the clustering signal of 
UHE protons. 
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suggestions and stimulating discussions. Anonymous referees are also acknowledged for 
their fruitful comments and suggestions. 

References 

[1] K. Greisen. End to the cosmic-ray spectrum? Phys. Rev. Lett., 16(17):748-750, 1966. 
[2] G.T. Zatsepin and V.A. Kuz'Min. Upper limit of the spectrum of cosmic rays. JETP Lett., 4, 
1966. 

[3] J. Abraham et al. Observation of the Suppression of the Flux of Cosmic Rays above 4xl0 19 eV. 

Phys. Rev. Lett., 101(6):61101, 2008. 
[4] R.U. Abbasi et al. First observation of the Greisen-Zatscpin-Kuzmin suppression. Phys. Rev. 

Lett, 100(10):101101, 2008. 
[5] J. Abraham et al. Measurement of the energy spectrum of cosmic rays above 1018 eV using the 

Pierre Auger Observatory. Phys. Lett. B, 685(4-5):239-246, 2010. 
[6] J.L. Puget, F.W. Stecker, J.H. Bredekamp. Photonuclear interactions of ultrahigh energy cosmic 

rays and their astrophysical consequences. Astroph. J., 205:638-654, 1976. 
[7] D. Allard, E. Parizot, A. Olinto, E. Khan, and S. Goriely. UHE nuclei propagation and the 

interpretation of the ankle in the cosmic-ray spectrum. Astron. Astroph., 443:L29~L32, 2005. 
[8] D. Harari, S. Mollerach, and E. Roulet. On the ultrahigh energy cosmic ray horizon. J. Cosm. 

Astrop. Phys., 2006:012, 2006. 
[9] D. Hooper, S. Sarkar, and A.M. Taylor. The intergalactic propagation of ultra-high energy cosmic 

ray nuclei. Astrop. Phys., 27(2-3) :199-212, 2007. 
[10] D. Allard and R.J. Protheroe. Interactions of UHE cosmic ray nuclei with radiation during 

acceleration: consequences for the spectrum and composition. Astron. Astroph., 502(3):803-815, 

2009. 

[11] D. Allard, N.G. Busca, G. Decerprit, A. Olinto, and E. Parizot. Implications of the cosmic ray 
spectrum for the mass composition at the highest energies. J. Cosm. Astrop. Phys., 2008:033, 
2008. 

[12] M. Kachelriefi, E. Parizot and D.V. Semikoz. The gzk horizon and constraints on the cosmic ray 
source spectrum from observations in the gzk regime. JETP Lett., 88(9):553-557, 2009. 

[13] M. Nagano and A. A. Watson. Observations and implications of the ultrahigh-energy cosmic rays. 
Rev. Mod. Phys., 72(3):689-732, 2000. 

[14] P. Bhattacharjee and G. Sigl. Origin and propagation of extremely high-energy cosmic rays. Phys. 
Rep., 327(3-4): 109-247, 2000. 



Influence of cosmological models on the GZK horizon of ultrahigh energy protons 23 



A.M. Hillas. The origin of ultra-high-energy cosmic rays. Ann. Rev. Astr. Astrop., 22(l):425-444, 
1984. 

C. T. Hill, D.N. Schramm, and T.P. Walker. Ultra-high-energy cosmic rays from superconducting 
cosmic strings. Phys. Rev. D, 36(4):1007-1016, 1987. 

V. Berezinsky and A. Vilenkin. Cosmic necklaces and ultrahigh energy cosmic rays. Phys. Rev. 

Lett, 79(26):5202-5205, 1997. 
V. Berezinsky, M. Kachehiefi, and A. Vilenkin. Ultrahigh Energy Cosmic Rays without Greisen- 

Zatscpin-Kuzmin Cutoff. Phys. Rev. Lett, 79(22) :4302-4305, 1997. 
A. Venkatesan, M.C. Miller, and A.V. Olinto. Constraints on the production of ultra-high-energy 

cosmic rays by isolated neutron stars. Ap. J., 484:323, 1997. 
G.R. Farrar and P.L. Biermann. Correlation between compact radio quasars and ultrahigh energy 

cosmic rays. Phys. Rev. Lett, 81(17):3579-3582, 1998. 

D. Fargion, B. Mele, and A. Salis. Ultra-High-Energy Neutrino Scattering onto Relic Light 
Neutrinos in the Galactic Halo as a Possible Source of the Highest Energy Extragalactic Cosmic 
Rays. Ap. J., 517:725, 1999. 

J. Arons. Magnetars in the metagalaxy: an origin for ultra-high-energy cosmic rays in the nearby 
universe. Ap. J., 589:871, 2003. 

E. Waxman, K.B. Fisher, and T. Piran. The Signature of a Correlation between Cosmic-Ray 
Sources above 10 19 eV and Large-Scale Structure. Ap. J., 483:1-7, 1997. 

G. Sigl, D.F. Torres, L.A. Anchordoqui, and G.E. Romero. Testing the correlation of ultrahigh 
energy cosmic rays with high redshift sources. Phys. Rev. D, 63(8):81302, 2001. 

Z. Fodor and S.D. Katz. Ultrahigh energy cosmic rays from compact sources. Phys. Rev. D, 
63(2):23002, 2000. 

P.G. Tinyakov and I.I. Tkachev. Correlation function of ultrahigh-energy cosmic rays favors point 

sources. JETP Lett, 74(l):l-5, 2001. 
J. Abraham et al. Correlation of the Highest-Energy Cosmic Rays with Nearby Extragalactic 

Objects. Science, 318(5852) :938, 2007. 
J. Abraham et al. Correlation of the highest-energy cosmic rays with the positions of nearby active 

galactic nuclei. Astrop. Phys., 29(3):188-204, 2008. 
P. Abrcu et al. Update on the correlation of the highest energy cosmic rays with nearby 

extragalactic matter. Astrop. Phys., 34(5):314-326, 2010. 
R.U. Abbasi et al. Search for correlations between HiRes stereo events and active galactic nuclei. 

Astrop. Phys., 30(4):175-179, 2008. 
A. Cuoco, R. D'Abrusco, G. Longo, G. Miclc, and PD Serpico. The footprint of large scale cosmic 

structure on the ultrahigh energy cosmic ray distribution. J. Cosm. Astrop. Phys., 2006:009, 

2006. 

E. Bugaev and P. Klimai. Interactions of high-energy cosmic rays with extragalactic infrared 

radiation background. Phys. Atom. Nucl, 70:156-158, 2007. 
J. Aublin et al. Discriminating models of UHECR sources with a log likelihood method. Nucl. 

Phys. B - Proc. Supp., 190:94-98, 2009. 
M. De Domenico et al. Bounds on the density of sources of ultra high energy cosmic rays from 

Pierre Auger Observatory data. Proc. 32nd LCRC, Beijing, 2011. 
D. Larson et al. Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 

Power Spectra and WMAP-Derived Parameters. Ap. J. S., 192(2):16, 2011. 
V. L. Ginzburg and S. I. Syrovatskii. The Origin of Cosmic Rays. Oxford: Pergamon Press, 1964. 
R. Aloisio and V. Berezinsky. Diffusive Propagation of Ultra-High-Energy Cosmic Rays and the 

Propagation Theorem. Astroph. J., 612:900, 2004. 
R. Aloisio and V.S. Berezinsky. Anti-GZK Effect in Ultra-High-Energy Cosmic Ray Diffusive 

Propagation. Astroph. J., 625:249, 2005. 
M. Lemoine. Extragalactic magnetic fields and the second knee in the cosmic-ray spectrum. Phys. 

Rev. D, 71(8):83007, 2005. 



Influence of cosmological models on the GZK horizon of ultrahigh energy protons 24 



R. Aloisio, V. Bcrczinsky, P. Blasi, A. Gazizov, S. Grigorieva, and B. Hnatyk. A dip in the 
UHECR spectrum and the transition from galactic to extragalactic cosmic rays. Astrop. Phys., 
27(1):76-91, 2007. 

V. Berezinsky and A.Z. Gazizov. Diffusion of Cosmic Rays in the Expanding Universe. II. Energy 
Spectra of Ultra-High Energy Cosmic Rays. Astroph. J., 669:684, 2007. 

R. Engel, D. Seckel, and T. Stanev. Neutrinos from propagation of ultrahigh energy protons. 
Phys. Rev. D, 64(9):93010, 2001. 

M. Ave et al. Cosmogenic neutrinos from ultra-high energy nuclei. Astrop. Phys., 23(l):19-29, 
2005. 

T. Stanev. High energy cosmic rays. Springer Verlag, 2009. 

G.R. Blumenthal. Energy loss of high-energy cosmic rays in pair-producing collisions with ambient 

photons. Phys. Rev. D, 1(6):1596-1602, 1970. 
J. P. Rachen. Interaction Processes and Statistical Properties of the Propagation of Cosmic Rays 

in Photon Backgrounds. PhD thesis, Bonn University, 1996. 
M.J. Chodorowski, A. A. Zdziarski and M. Sikora. Reaction rate and energy-loss rate for photopair 

production by relativistic nuclei. Astroph. J., 400:181-185, 1992. 

F. W. Stecker. Effect of photomeson production by the universal radiation field on high-energy 
cosmic rays. Phys. Rev. Lett., 21(14):1016-1018, 1968. 

J. P. Rachen and P.L. Biermann. Extragalactic Ultra-High Energy Cosmic-Rays-Part One- 
Contribution from Hot Spots in Fr-II Radio Galaxies. Astron. Astroph., 272:161, 1993. 

L.A. Anchordoqui, M.T. Dova, L.N. Epele and J.D. Swain. Effect of the 3 K background radiation 
on ultrahigh energy cosmic rays. Phys. Rev. D, 55(12):7356-7360, 1997. 

T. Stanev. Propagation of ultrahigh-energy cosmic rays. New J. Phys., 11:065013, 2009. 

J. P. Vallee. Cosmic magnetic fields as observed in the universe, in galactic dynamos, and in the 
milky way. N. Astron. Rev., 48(10):763-841, 2004. 

TR Seshadri and K. Subramanian. Cosmic microwave background bispectrum from primordial 
magnetic fields on large angular scales. Phys. Rev. Lett., 103(8):81303, 2009. 

P. Trivcdi, K. Subramanian, and TR Seshadri. Primordial magnetic field limits from 
cosmic microwave background bispectrum of magnetic passive scalar modes. Phys. Rev. D, 
82(12):123006, 2010. 

D. Ryu, H. Kang, J. Cho, and S. Das. Turbulence and magnetic fields in the large-scale structure 
of the universe. Science, 320(5878):909, 2008. 

S. Das, H. Kang, D. Ryu, and J. Cho. Propagation of ultra-high-energy protons through the 

magnetized cosmic web. Astroph. J., 682:29, 2008. 
J.E. Everett, E.G. Zweibel, R.A. Benjamin, D. McCammon, L. Rocks, and J.S. Gallagher III. The 

milky way's kiloparsec-scale wind: A hybrid cosmic-ray and thermally driven outflow. Astroph. 

J., 674:258, 2008. 

J.E. Everett, Q.G. Schiller, and E.G. Zweibel. Synchrotron constraints on a hybrid cosmic-ray 
and thermally driven galactic wind. Astroph. J., 711:13, 2010. 

E. Armengaud, G. Sigl, T. Beau, and F. Miniati. CRPropa: A numerical tool for the propagation 
of UHE cosmic rays,[gamma]-rays and neutrinos. Astrop. Phys., 28(4-5):463-471, 2007. 

A. Neronov and I. Vovk. Evidence for strong extragalactic magnetic fields from fcrmi observations 

of tev blazars. Science, 328(5974):73, 2010. 
M. Ahlers and J. Salvado. Cosmogenic gamma rays and the composition of cosmic rays. Phys. 

Rev. D, 84(8):085019, 2011. 

G. Decerprit and D. Allard. Constraints on the origin of ultra-high-energy cosmic rays 
from cosmogenic neutrinos and photons. Astron. Astroph., 535:A66, 2011. 

V. Berezinsky, A. Gazizov, and S. Grigorieva. On astrophysical solution to ultrahigh energy cosmic 

rays. Phys. Rev. D, 74(4):043005, 2006. 
W. H. Baumgartner, J. Tueller, C. Markwardt, and G. Skinner. The Swift-BAT 58 Month 

Survey. In Bulletin of the American Astronomical Society, volume 42 of Bulletin of the American 



Influence of cosmological models on the GZK horizon of ultrahigh energy protons 25 

Astronomical Society, page 675, 2010. 
[65] M. De Domenico, A. Insolia, H. Lyberis, M. Scuderi. Multiscale autocorrelation function: a new 

approach to anisotropy studies. J. Cosm. Astrop. Phys., 2011:008, 2011. 
[66] S. Kullback and RA Leibler. On information and sufficiency. The Annals of Mathematical 

Statistics, 22:79-86, 1951. 



